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Abstract 

The relation between rainfall and water accumulated in reservoirs 
comprises nonlinear feedbacks. Here we show that they may generate 
alternative equilibrium regimes, one of high water-volume, the other 
of low water-volume. Reservoirs can be seen as socio-environmental 
systems at risk of regime shifts, characteristic of tipping point transi¬ 
tions. We analyze data from stored water, rainfall, and water inflow 
and outflow in the main reservoir serving the metropolitan area of Sao 
Paulo, Brazil, by means of indicators of critical regime shifts, and find 
a strong signal of a transition. We furthermore build a mathematical 
model that gives a mechanistic view of the dynamics and demonstrates 
that alternative stable states are an expected property of water reser¬ 
voirs. We also build a stochastic version of this model that fits well 
to the data. These results highlight the broader aspect that reservoir 
management must account for their intrinsic bistability, and should 
benefit from dynamical systems theory. Our case study illustrates the 
catastrophic consequences of failing to do so. 
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1 Introduction 


Complex socio-ecological and socio-environmental systems often involve in¬ 
teractions between natnral elements and tinman action in a nonlinear and 


adaptive way (! 

^olke, 

2006 

). Such systems have received much attention 

(Gordon et al. 

2008 

|Osterblom et al. 

2013 

) recently and may display 


many of the common dynamical featnres present in natnral systems. Here 
we will be interested in the existence of alternative stable states, representing 
different possible dynamical regimes, and transitions between them. These 
tipping-point transitions have been well stndied and are common in natnral 


systems ( 

Scheffer et al. 

2001|. E 

et al. 

2004 

Scheffer & Carpenter 

& van Nes, 

2007 

Carpenter 

2005 


Desertihcation (Klansmeier, 1999 


2003) and lake entrophication (Scheffer 


Meron 


a plethora of cases. Characterization of a transition from time-series data 
can be assessed nsing indexes related to either retnrn-time to eqnilibria or 
variability near the tipping point (Dakos et al. , 2012). These techniqnes. 


together with models that encompass the main featnres of the system, may 
be extended to systems where the interaction with the hnman factor creates 
the possibility of a certain level of control or, at least, of action attempting 
to bring the system to a sitnation that is considered desirable. Onr atten¬ 
tion will be focnsed on the dynamics of water reservoirs, which has not been 
stndied nnder this light, and on the characterization of alternative stable 
states and associated transitions. We show a case in which these dynamic 
properties had been neglected in the management of one of the largest water 
reservoirs in the world, the Cantareira system, which serves the Metropolitan 
Area of Sao Panlo (MASP). 

MASP faces an nnprecedented crisis of water snpply. This exceptional 
sitnation can be readily seen from the plot of the volnme of water stored 
the Cantareira system, the main reservoir serving the area, over the past 
twelve years (Fig.([^). The volnme of water decreased sharply from mid- 

2013 and the operational capacity of the reservoir was depleted in July 2014. 
Since then water withdrawal is done by pumping of the so-called “strategic 
reserve” or “dead volume”. The Sao Paulo Water Company, SABESP, began 
to reduce withdrawals in January 2014, and by May 2015 the total outflow 
was 40% of the average values. Additionally, the last rainy season (October 

2014 - March 2015) provided more rainfall compared to the previous two 
years. Despite that, only of 15% of total volume had been recovered, and 
the reservoir remains operationally exhausted. This situation seems typical 
of catastrophic regime shifts, driven by a tipping-point transition. In this 
work, we will establish this regime transition on hrm grounds and discuss its 
consequences. For sake of brevity, we will use the term reservoir to mean a 
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system of interconnected reservoirs. 
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Figure 1. Rainfall, water flows and volume stored in the Cantareira 
system since 2004. Upper panel: notwithstanding the inter-annual trend, a 
clear seasonal fluctuation is present in the rainfall (gray) which was followed by 
the volume of stored water (percentage of operational volume, blue line) until 
mid 2013. From this point on, the volume of water decreased to very low levels, 
without any seasonal trend. A recovery of small amplitude is seen from Feb- 
May 2015 due to reduction in outflow. This is shown in the lower panel: the 
water inflow into the reservoir (gray) showed seasonal fluctuations coupled to 
rainfall. The reservoir operators kept outflow (blue line) around the maximum 
average allowed (36 exceeded only to avoid overflow). Operators started 

a gradual decrease in outflow in the beginning of 2014, when the operational 
volume was close to depletion. Since the depletion (July 2014), water is available 
only by pumping, and outflow reduction proceeded in a slow pace till January 
2015, when it was abruptly cut by 25%. 


We propose that reservoir water-volume is subject to drastic regime changes 
due to the underlying bistability of the system. Bistability is the condition 
in which there are at least two possible stationary states of a system for the 
same set of parameters. In our case, the parameter is rainfall and the vari¬ 
able that may assume two possible alternative stable states is the volume of 
water stored in the reservoir. If the reservoir is at high levels the catchment 
has more stored water and thus more of the rainfall will flow into water the 
reservoir. On the other hand, for low water levels the most obvious change 
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is that much of the water is absorbed by dry soil and thus a lower propor¬ 
tion of rainfall becomes stored water in the reservoir. Many other factors 
may be play and the exact relation between rainfall and the amount of water 
stored in the reservoir is an instance of the rainfall-runoff problem in hydrol¬ 
ogy (Beven, 2011). What matters for our discussion is that these feedback 
mechanisms are drivers of critical transitions. 

Both states, high-level and low-level regime, are resilient in ranges of rain¬ 
fall and outflow parameters that may overlap, leading to hysteretic behavior: 
the paths to and from the low-level state are not the same. Resilience of the 
high-level water volume regime may be lost, for instance, if the rainfall index 
falls below a certain critical value, or if the reservoir is over-exploited. The 
system then transits to the low-level regime, in a tipping-point transition 
entailing a regime shift. The dramatic side of this fact is that the low-level 
regime is also resilient, so either rainfall index has to increase beyond the 
mentioned critical value, or outflow has to be reduced drastically in order to 
produce a backwards transition to the high-level regime. This state of af¬ 
fairs prompts us to characterize transitions by indexes related to the critical 


slowing down phenomenon (Dakos et al. , 2015). These indexes are usually 


connected to research of early-warning signals of critical transitions. Discus¬ 
sions in recent literature have raised doubts about how early the signals are 
dBoettiger et al. 2013), but this is not actually relevant in the present case. 


as we are concerned mainly with characterization and not anticipation of the 
transition. 

In what follows we proceed by the following steps: we hrst inspect the 
data to establish the existence of a feedback mechanism; next we look for 
indicators of critical transitions in two steps, hrst by a qualitative assessment 
of the data, and second by demonstrating the existence of a peak in the vari¬ 
ance characteristic of near tipping point regions. We then turn to a further 
step by proposing a mathematical model for the dynamics of the system, 
in two versions. First, a deterministic version where we can understand the 
dynamical behavior of reservoirs. Next, we propose an stochastic extension 
of the model, via a stochastic differential equation which can be htted to 
actual data accurately and used for short-term extrapolations. 


2 Model Independent Results 

The feedback mechanism 

Given a certain amount of rainfall, how much water will effectively how into 
the reservoir? This is connected to the rainfall-runoh problem, central in 
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hydrology. The big picture is that the wetter the catchment the more of 


the rainfall is converted in runoff water (Beven, 2011), that in this case 
will flow to the reservoir. A complete assessment of this situation involves 
modeling the rainfall-runoff relation in an explicit spatial setting, considering 
soil heterogeneities, terrain and many other variables. However, in order to 
proceed in a simpler way, we postulate that the volume of water stored in 
the reservoir is a surrogate for the drainage basin condition. To verify that 
this is a sound assumption, we show in Fig.(|^ the ratio of water inflow to 
the rainfall in the Cantareira reservoir as a function of the volume, yielding a 
clear trend: the higher the volume, the more efficiently the rainfall becomes 
water flowing into the reservoir. 



5 10 25 50 100 

Stored water (% maximum volume) 

Figure 2. Dependence of the ratio inflow to rainfall on the stored wa¬ 
ter in the reservoir. The Cantareira reservoir shows the well-known non-linear 
transfer of rainfall to runoff. The volume of water in the reservoir is a surrogate 
of total water in the whole catchment. The higher the volume the wetter the 
catchment and the larger the inflow/rainfall ratio. The data is the same of Fig. 

El 
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Assessment of a critical transition 

In Fig.0 we show the data for the Cantareira reservoir in the rainfall vs. 
stored-water plane. In normal sitnations, we wonld expect to have nearly 
closed curves in the upper part of the plane. These curves come from the 
seasonality of rainfall, which induces similar oscillations in the water stored 
in the reservoir. However, what is seen in the plot is that from 2013 on, the 
system spins down to a new cycle on the lower part of the plane. This lower 
cycle represents, again, seasonal variations, but the volume now oscillates 
around much smaller values. In the inset of Fig.([^ we show what would be 
the result obtained from an ordinary differential equation displaying bistable 
behavior, in which the rainfall parameter oscillates in time and has a decreas¬ 
ing trend, inducing a regime change. Similarities are evident. The specihc 
equations that generate Fig.([^ will be discussed in the Mathematical Models 
section. 



Figure 3. Observed and simulated transition between two dynamical 
regimes in the rainfall-volume plane. In the main panel we show the actual 
data for the Cantareira system. In the inset we show the results of a differential 
equation describing this system, which exhibits bistable behavior. The equation 
was simulated with variation of rainfall with a time oscillating part plus a long¬ 
term decreasing trend. 
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The previous discussion gives a qualitative view of a transition, but we 
would like to take a step forward and provide a more quantitative assessment. 
In recent years, a toolbox for detecting early-warning signals of critical tran¬ 
sitions has been developed (Dakos et al. , 2012 Dakos, 2015). These are 


based on the critical slowing down phenomenon, implying that the lag-auto- 
correlation function and variance indicators should increase near the transi¬ 
tion. As discussed in Boettiger & Hastings (2012), a series of underpinnings 
exist concerning the reliability of these indicators, which may fail to antici¬ 
pate the occurrence of a transition. However, in this work we are interested 
in characterizing the regime shift rather than anticipating it. In order to do 
so, we evaluated evidence of a critical transition with a diffusion-drift-jump 
model {ddj). If the time series is long enough and has high frequency of ob¬ 
servations, this model gives an accurate description of the underlying process 
that generates the data (Carpenter & Brock, 2011). Variations along time 
are described as a combination of instantaneous changes by deterministic 
trends (the drift), instantaneous variation (diffusion) and occasional uncor¬ 
related changes (jumps). The model also allows to estimate the variance 
conditioned to these components. In theory, this conditional variance goes 


to inhnity at bifurcation points of the underlying dynamics (Dakos et al. 


2012). We found a very clear peak in the conditional variance near the be¬ 
ginning of 2014 (Fig.(|^). This gives the hrst characterization of an indicator 
of critical regime shifts for water reservoirs. Running-window statistics like 


2012) followed the 


variance and autocorrelation at lag one ( |Dakos et al. 
same trend. 

As discussed in Dakos et al. (2015), critical transitions are better assessed 
when, together with a statistical indicator, like the conditional variance used 
in Fig.(|^, a mechanism leading to bistability is evident. The feedback dis¬ 
cussed above is the basic mechanism we hypothesize. A further step, there¬ 
fore, is to build a mathematical model that takes into account the feedback 
and may be fitted to the data, and show that it exhibits bistability. This is 
what we present in the next section. 


3 Mathematical models for reservoir dynam¬ 
ics 

Deterministic model 

The change in time of stored-water volume in a reservoir is determined by 
inflow and outflow rates. Outflow rates are a consequence of management 
decisions. Inflow comes ultimately from rainfall through a set of hydrological 
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Figure 4. Assessment of critical regime shift in a water reservoir 
by a peak in the conditional variance. Conditional variance from a drift- 
diffusion-jump model (ddj) for the Cantareira water volume as a function of rel¬ 
ative stored volume (left), and of time (right). The ddj model was fitted to 
the daily time series of the relative stored volume, log-transformed as recom¬ 


mended by Dakos et al. (2012). The time series (right panel) shows an abrupt 


increase of ddj conditional variance at the beginning of 2014, which indicates 
that a state transition occurred. The left panel shows that the high variance 
values occurred in a narrow range of low volumes in the reservoir, recorded in 
the end of 2013 and beginning of 2014. 


processes. This leads to an equation of the form: 

'^ = -s[V)+r{V,R) , (1) 

where V is the volume of water stored in the reservoir, s(F) is the outflow 
rate and r(l/, R) is the inflow rate, with R a measure of the rainfall. We now 
discuss these terms separately. 

Outflow. We assume that withdrawal of water is managed and saturates 
to a maximum, sq; when the stored volume is high, but is reduced as water 
becomes scarce (which could be done both to prevent complete disaster or 
due to the difficulty in pumping water from low levels). Such assumption is 
realistic for the Cantareira system: while stored volume was above 33% water 
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withdrawals oscillated around 36 according to the maximum quota set 

by the Brazilian Water Agency, ANA (Fi^^. As the volume dropped outflow 
decreased because ANA pushed the main reservoir operator SABESP to cut 
down withdrawals. To model this situation we use the following simple form: 


s(E) 


K + V ' 


( 2 ) 


where K is the volume of water in which the withdrawal is reduced to so/2. 
It is a measure of how cautious the water management is. 


w 



Q 





0 40 80 

Stored water (% operational volume) 


Figure 5. Water withdrawal from the Cantareira reservoir depends 
on stored water volume. Withdrawals by the reservoir operators saturated to 
a maximum average value when stored water was above one third of operational 
capacity. Note also that when stored water fell below -15% operators plunged 
withdrawals from ca. 23 to 14 vN‘s~^. Since then releases are kept at these lev¬ 
els despite an increase in the stored volume, to save the current volume (darker 
points). 
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Inflow. As already implied in the notation r(l/, R), inflow depends not only 
on rainfall, but also on volume. As discussed above, we propose that inflow 
is higher for higher volumes. A simple assumption is thus: 


r{V,R) 


aRf^V^ ,ifV<Vma. 
0 , otherwise 


(3) 


where a, 0 < /3 < 1 and 0 < 7 < 1 are phenomenological parameters. 
We also introduced a maximum volume for the reservoir, Vmax- Once it is 
attained, no increase in stored water is possible and excess inflow is released 
downstream through a spillway. 


Putting everything together, we come to the following differential equation: 


dV 

dt 


soV 
K + V 


+ aRf^V^ , 


(4) 


if P < Vmax- If P > Vmax, the secoud term on the right-hand side is zero. 
Equation (|^ is not meant to be a full model of all hydrological processes, 
but rather a description of the basic dynamical factors affecting the volume 
of water stored in a reservoir. In spite of its simplicity, a stochastic version 
of it fits very well to observed data, as will be seen below. 

If the rainfall is constant over time. Equation (|^ has either one or two 
stable fixed points, depending on the value of R. Fig.(|^ shows the general 
situation, assuming specific values for the parameters. The existence of a 
bistability region is clear. 

The existence of the bistability region depends essentially on the fact 
that inflow is dependent on the volume. The form of the outflow term is not 
essential: all that is required is that it goes to zero as P goes to 0 and that 
it saturates at some value. The existence of a maximum volume is, from the 
mathematical viewpoint, a form of avoiding arbitrarily large volumes. 

Obviously, in real situations R is not constant. It has clear seasonal 
variations and may be subject to long-term trends. Solutions, in this case, 
will not tend to a constant value. Bistability will be reflected in the ex¬ 
istence of different regimes, corresponding to oscillations around the fixed 
points of the time-independent problem, as long as the seasonal variations 
and trends are not too large. In Fig.([^ an example is shown, with R = 
Ro{l + asm(27ift)){1 — bt). We also plot the solution as a curve in the VxR- 
plane (inset of Fig. ([^), showing the distinctive feature of oscillations in the 
upper part, spinning down on a short time scale to the lower part. This pat¬ 
tern is present also in the observed relationship between rainfall and stored 
volume in the reservoir (main panel of Fig.([^. 
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Figure 6. Fixed points of Equation Q as a function of rainfall R. We 

assumed that Vmax = 1; and V is measured in fractions of the maximum volume. 
Parameter values are: a = 1, so = 1, R = 0.2, /3 = 1, 7 = 2/3. The region of 
bistability is shadowed. 


Stochastic model 

In order to compare theory to data we need to take a further step and allow 
for stochastic fluctuations. This leads us to consider the following stochastic 
differential equation: 

dV = {aR^ - E)dt + aVdW , (5) 

where E is the outflow of water R is the mean rainfall in the previous 

30 days (mm) and W is the standard Wiener process. The term aVdW 
expresses an instantaneous stochastic Gaussian noise in the stored volume, 
with zero mean and standard deviation proportional to the stored volume. 
We also assume that the observed values of stored volume follow a Gaussian 
distribution with mean value equal to the expected value of the stochastic 
equation and an unknown standard deviation e. The reason to take R as 
a 30-day mean is empirical: volume does not respond instantaneously to 
rainfall, and complex hydrological processes tend to spread effects over time. 
This non-locality in time is encompassed by taking a mean rainfall over a 
certain period. The period itself (30 days) is justihed a posteriori as the one 
giving the best £t to data. 

Equation ([^ keeps the deterministic skeleton from the previous discus¬ 
sion, at the same time allowing us to take V, R, E and inflow from time-series 
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Figure 7. Regime shift in our mathematical model. The volume time 
series is a representative numerical solution of equation For this solution rain¬ 
fall is assumed as R{t) = -Ro(l + asin(27r/t))(l — bt), with Rq = 1.1, a = 0.4, 

/ = 0.4 and b = 0.01. All the other parameters are the same as in Fig.®. The 
parametric plot V{t)K.R{t) for this solution is given in the inset of Fig.([3F 


of recorded data for the Cantareira reservoir. We recall that our model pre¬ 
dicts catastrophic shifts because water inflow depends on the stored volume, 
which creates a feedback mechanism. To evaluate the support provided by 
data to this hypothesis, we evaluated a competing model in which inflow 
depends only on rainfall: 

dV = {aR^ - E)dt + aVdW ( 6 ) 

The dataset are the time series of daily inflow, outflow, rainfall, and stored 
water volume. The model provided good hts for time spans up to one year. 
In this paper we focused on the last 365 of the time series (May 31 2014 to 
May 31 2015). Because simultaneous htting of all hve parameters ( a, /?, 7 , 
a and e) did not converge or provided unreasonable estimates, we callibrated 
both models in 3 steps: 


1. The observation error e was conservatively estimated from trajectory 


matching (King et al. , 2009) 


2. The exponent 7 was estimated as the slope of a Gaussian linear regres¬ 
sion of reservoir area in function of reservoir volume in log scale. 
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3. Parameters a, f3 and a were estimated with Bayesian particle filter (Lin 
|fc West 2001) modified by King et al. (2009). 


Equation!^ which describes expected water inflow as a fnnction of rainfall 
and current stored volume, provided a much more plausible £t to the observed 
time series (log-likelihood ratio: 12.11, Tab. [^. The model that does not 
take into account the effect of volume on inflow (Eq.([^) underestimated the 
stored volume in most of the period, and did not predict the increase and 
further stabilization of stored volume since February 2015 (Fig.(|^). The 
better £t of equation supports the hypothesis that the ratio of inflow to 
rainfall depends on the volume. This feedback is caused by the interaction 
between rainfall and stored volume, a surrogate of the hydrological state of 
the catchment. This, in turn, substantiates our statement about the existence 
of alternative states due to a feedback process. 


Table 1. Estimates of parameters and Log-likelihood (LL) for two competing 
models to predict the changes in stored water volume in Cantareira Reservoir. 
The model that describes water inflow as a function of rainfall and stored water 
(/(rain, volume), Eq.([^) had a much higher likelihood and thus provides a much 
more plausible description of the time series. 


Model 

a 


7 

a 

e 

LL 

Eq^ 

) /(rain, volume) 

5.998 

1.043 

0.590 

0.00231 

3.2 X 10^ 

-6676.0 

Eq.i 

5 /(rain) 

392494.3 

0.926 

- 

0.00948 

3.2 X 10^ 

-6688.1 


4 Concluding remarks 

We have demonstrated that the dynamics of the volume of water stored in 
a reservoir is subject to abrupt regime shifts and hysteretic behavior. These 
results follow from: (i) a qualitative assessment of the dynamics; (ii) a quan¬ 
titative evaluation of indexes characterizing regime shifts; (iii) a modeling 
approach that gives a mechanistic view of the dynamics; and (iv) a £t of a 
corresponding stochastic model to the data. 

These results have serious signihcance for water reservoirs management. 
Resilience of the normal operating conditions may be lost under drought 
conditions in favor of a new low water levels regime, which is itself resilient 
under a range of external conditions. Recovery of the normal condition can, 
thus, be difficult, and may demand potentially extreme measures. These 
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Figure 8. Time series of observed and predicted stored water volume 
in the Cantareira reservoir system. The black line shows the observed time 
series which are daily records from 2014-05-31 to 2015-05-31. The orange lines 
and shaded areas show the expected trajectories and Bayesian credibility interval 
from a mathematical model of observed inflow as function of rainfall (Equation 
1^. The blue lines and shaded areas show the expected values and intervals from 
a model of inflow as function of both rainfall and the stored volume (Equation 
[^. Trajectories and intervals were calculated from 5000 numerical simulations 
of each model, with parameters sampled from the posterior distribution of the 
Bayesian particle filter estimation. 


measures are those that allow for decrease in the outflow of the reservoir, 
representing a possible burden to populations served by the reservoir. This 
implies that cautious management of a reservoir should avoid the regime 
shift. The major question then becomes: how to achieve such management? 
We have taken a first step here, identifying the need to acknowledge the 
potential for such catastrophic transitions. Under this perspective, early 
warning signals of catastrophic regime shifts are a possible avenue to tackle 
the problem in the absence of detailed hydrological data, but the extent of 
their predictive power has yet to be ascertained. 

Studies in water resources have mainly favored full-blown models to depict 
the ever increasing understanding of hydrological processes ( Beven| 2011). 
Nevertheless simple mathematical descriptions of the hrst principles that 
operate in catchments also proved to be useful both in fitting data and un- 
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derstanding the underlying mechanisms (|Beven, 2011 Young et al. , 2004). 


We followed this synthetic approach to propose a model akin to those used 
in system dynamics. Our models are a bare-to-the-bones description on how 
rainfall turns into stored water in a reservoir. In doing that we depict many 
hydrological processes in a few hxed parameters. What could remain from 
such a phenomenological approach are very general dynamic properties. In 
this sense, we show how alternative states and regime shifts result from a very 
simple hrst principle, namely, nonlinear feedback mechanisms. We show that 
a model that has these feedbacks - and not much more than that - provides 
a more plausible description of observed data. In doing so we aim to high¬ 
light key properties that water reservoir share with other socio-environmental 
systems, which can contribute to fully-fledged models. 

In this study we focused on a specihc reservoir, the Cantareira system 
in the metropolitan area of Sao Paulo. Let us now elaborate on this spe¬ 
cihc case. Signals of below the average rainfall were present since mid-2013, 
with corresponding decrease in water inhow. The conditional variance in¬ 
creased steadily in the very beginning of 2014, indicating an ongoing regime 
shift. The average inhow rate in January 2014 was 15.7 m^/s, whereas the 
outhow was kept at 34.3 m^/s, a very unusual dehcit for the rainy season 
in the region. Actually, the 30-day moving average of inhows and outhows 
presented dehcits continually from May 2013 to February 2015, contrary 
to the recurrent pattern where the period from December-May always had 
greater inhows than outhows. The operators began decrementing releases 
only in March 2014, and at a pace that was not enough to bring back the 
system to the normal operating conditions. The obvious question is why the 
management policies could not avoid the collapse of the system. Official doc¬ 
uments available to the publi c ( Age nda Nacional de Aguas & Departamento 


de Aguas e Energia Eletrica 


2015) show that the Brazilian Water Agency 


(ANA) and the Sao Paulo State Water Department (DAEE) agreed that op¬ 
erators {SABESP and other water companies) could take up to 36 m^s—1 of 
water from the Cantareira system. The agreement also included a reservoir 
rule curve that prescribed maximum allowed withdrawals according to the 
stored volume of the system. These were not respected since January 2014. 
Even if the limits were met, they were probably too loose. For instance, the 
outflow allowed for the rainy months of November/December 2014 would be 
27 m^s—1 even if the operational capacity had been completely depleted. 
Since March 2014, ANA and DAEE abandoned the rule curve and adopted 
tighter limits. The maximum values have been negotiated periodically, and 
are currently 17.5 m^s—1 (June-August 2015) and 13.5 m^s—1 (September- 
November 2015). Nevertheless, the National Center for Surveillance of Natu¬ 
ral Disasters (http://www.cemaden.gov.br/) forecasts that the Cantareira 
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system will stay below its operational capacity at least until the beginning 
of 2016. 

In summary, our results and our study case show that the management 
of reservoirs should take alternative regimes into account and avoid a tran¬ 
sition to low-volume regimes. Failing to do so represents a prolongated bur¬ 
den, extending well beyond the period of anomalous rainfall, because out¬ 
flow has to be kept as low as possible until a backwards transition occurs. 
Therefore managers should act as another feedback mechanism in the socio- 
environmental system that keeps it in the desired regime notwithstanding 
external forces like climate anomalies. In not doing so, managers of the 
Cantareira system acted like one more external force that pushed the reser¬ 
voir to a catastrophic shift. 


Data Sets 

Data for the time series of stored water volume in the main reservoir system 
of Sao Paulo State (Cantareira) and for the corresponding rainfall index were 
obtained from the Sao Paulo Water Company (SABESP) database, which can 
de accessed at http://www2.sabesp.com.br/mananciais/DivulgacaoSiteSabesp. 
aspx. Data for the inflow and outflow were obtained from http://www2. 
sabesp.com.br/mananciais/Divulgacaopcj.aspx and processed through 
Optical Character Recognition (OCR) software. Starting January 15th, 2015, 
a more convenient tool has been provided by SABESP at http://site, 
sabesp.com.br/site/interna/Default.aspx?secaoId=553, All data and 
codes used in the analyses are available at https: //github. com/cantareira/ 
plos. 
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